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SECTION  I 
INTRODUCTION 


The  aim  of  this  study  has  been  to  compare  the  expected  perform¬ 
ance  of  several  instrumentation  systems  for  measuring  transmission  at 
94  GHz  under  various  weather  conditions.  This  frequency  was  selected 
because  it  is  one  of  the  “windows"  of  low  molecular  absorption  as  shown 
in  Figure  1.  The  problem  of  calculating  the  attenuation  through  rain 
and  fog  due  to  scattering  and  absorption  in  terms  of  meteorologically 
measurable  quantities  was  also  attacked  as  a  part  of  the  system  evalu¬ 
ation  problem. 

The  following  receiver  systems  were  studied: 

1.  Square-Law  Detection, 

2.  Quadrature  Receiver, 

3.  Post-Detection  Integration. 

The  receiver  with  a  square-law  detector  was  considered  both  in  a  one¬ 
way  transmission  system  and  in  a  two-way  transmission  system  with  the 
use  of  a  calibrated  target,  i.e.,  as  a  radar  transmissometer .  In  both 
cases,  the  receiver  signal  was  considered  to  be  chopped  at  an  audio 
frequency.  The  quadrature  receiver  was  studied  for  the  reception  of 
both  continous-wave  and  pulsed  signals  in  the  radar  transmissometer 
configuration.  The  receiver  which  employs  post-detection  integration 
was  considered  to  be  used  with  a  pulsed  transmitter:  a  standard  pulse 
radar  configuration. 

The  signal -to-noise  ratio  expressions  and  an  expression  for 
attenuation  by  rain  and  fog  are  given  in  Section  H.  The  results 
obtained  from  computer  programs  utilizing  these  equations  are  contained 
in  Section  III. 
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Figure  1.  Specific  attenuation  in  clear  air  versus 
frequency  (after  Rosenblum  [1]). 
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SECTION  II 

THE  SYSTEM  PERFORMANCE  EQUATIONS 


The  range  at  which  an  output  signal-to-noise  ratio  of  10  dB  would 
be  realized  was  taken  as  a  criterion  of  system  performance.  Therefore, 
the  aim  of  this  section  is  to  derive  signal-to-noise  ratio  expressions 
for  the  five  types  of  instrumentation  which  were  presented  in  Section  I. 
As  a  first  step,  the  signal-to-noise  ratio  expressions  for  the  receiver 
systems  will  be  derived  at  a  point  before  any  detection  or  non-linear 
signal  processing  is  employed,  i.e.,  at  the  output  of  the  inter¬ 
mediate  frequency  amplifier.  These  expressions  depend  on  whether  a 
one-way  path  or  a  two-way  (radar)  mode  is  considered,  but  they  do  not 
depend  on  the  type  of  signal  detection  employed. 

A.  Derivation  of  Pre-Detection  Signal - 

to-Noise  Ratio  Expressions 

Each  system  falls  into  one  of  the  following  two  groups: 

1.  The  system  utilizes  the  transmission  of  a  signal 
over  a  one-way  path,  or 

2.  the  receiver  is  located  physically  close  to  the 
transmitter  and  the  signal  from  the  transmitter 
is  reflected  from  a  calibrated  target  to  the 
receiver.  In  essence,  this  is  a  radar  system 
although  without  a  requirement  for  range  resolution. 

For  systems  of  the  first  kind,  the  power  density  at  the  receiver 
antenna  is  given  by 
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* 


(1) 


ptGt 


4tR 


where  is  the  transmitted  power,  G^  is  the  transmitter  antenna  gain 
in  the  direction  of  the  receiver,  and  R  is  the  distance  between  the 
transmitter  and  receiver  antennas,  all  in  consistent  units.  Then,  the 
received  power  is  given  by 


P 

r 


4-R2  e 


(2) 


where  Pr  is  the  power  received  by  the  receiving  antenna  and  A  ,  in 
units  of  area,  is  the  effective  aperture  of  the  receiving  antenna, 
given  by 


A 


e 


(3) 


where  is  the  gain  of  the  receiver  antenna  in  the  direction  of  the 
r 

transmitter,  and  X  is  the  wavelength  of  the  signal.  Thus,  from 
Equations  (2)  and  (3),  one  obtains  the  transmission  equation  given 

by  [2] 


p 

r 


ptGtV 

(4*)2R2 


(4) 


If  there  are  attenuating  factors  over  the  path  of  propagation 
such  as  rain  or  fog,  and  if  there  is  loss  in  the  transmitter  and 
receiver  systems,  then  the  inclusion  of  these  effects  leads  to 


P  -  PtGtGrx  Vr 

r  (4tt)2R2 


,-ocR/l  0 


(5) 


where  a  is  the  specific  attenuation  in  decibels  per  unit  length,  and 
lt  and  Lr  are  transmitter-to-antenna  and  antenna-to-receiver  component 
losses,  respectively. 
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The  thermal  noise  power  at  the  output  of  the  pre-detection 
filter  (usually  simply  the  intermediate  frequency  amplifier)  is 


.  =  kToBF  ,  (6) 

-23 

where  k  is  Boltzmann's  constant  (1.38*10  joules/K),  B  is  the 

effective  receiver  pre-detect ion  filter  bandwidth  (Hertz),  T  is  the 

o 

reference  noise  temperature  (290  K).  and  F  is  the  noise  figure  of  the 
receiver  referred  to  290  K. 


Thus,  from  Equations  (5)  and  (6),  one  obtains 

(S/N)  =  tG^r?--L-—  io"‘R/1°  .  ( 

(4-r )  FrkTQBF 

This  is  the  s ignal -to-noi se  ratio  at  the  output  of  the  intermediate- 
frequency  amplifier  for  one-way  transmission. 


For  systems  of  the  radar  type,  the  power  density  at  the  target 
is  given  by 

p(  .  Wi  .  (! 

r  4t(T 


The  power  reradiated  in  the  receiver  direction  is 

p.  ,  Wtl,0-«R/io 

4^R^ 


(9) 


where  a,  in  units  of  area,  is  the  radar  cross  section  of  the  target. 
The  power  density  of  the  echo  signal  at  the  receiver  antenna  is  given 
by 


P  = 


PtGtLta  i n-2aR/l 0 
?  ?  u 
(4ttR  *-)L 


00) 


and  the  received  power  is  obtained  by  multiplying  p  by  the  effective 
area  of  the  aperture,  leading,  by  use  of  Equation  (3),  to  the  radar 
equation  including  losses  [2], 


PtGt.V  jLtLr  1  -2.R/10 
(4-)V 


(11 


Thus,  using  Equations  (6)  and  (11),  the  pre-detection  signal-to- 
noise  ratio  for  the  radar  transmi ssometer  is 


,r;ilX  PtGtV  lLtLr  ,  n-2  xR/1 0 
(S/N)  =  - -  10 

(4- )  Jf<  kT  BE 
o 


(12) 


Equations  (7)  and  (12)  give  the  pre-detection  s ignal -to-noi se 
ratio  expressions  for  the  two  transmitter-receiver  geometries. 

B.  Processed  Signal -to-Noise  Ratio 
Expressions 

In  this  section,  the  five  systems  introduced  in  Section  I  will  be 
discussed,  and  for  each  one  the  excression  for  signal -to-noise  ratio, 
after  detection  and  signal  processing  has  been  employed,  will  be  given. 


1  .  One-way  transmission  and 
square-law  detection 

This  system  transmits  and  receives  a  chopped  signal  over  a  one¬ 
way  path.  At  the  receiver,  the  signal  is  heterodyned  to  an  intermedi¬ 
ate  frequency,  amplified,  and  detected  by  a  square-law  device  which  is 
followed  by  a  narrow-band  (audio)  filter  as  shown  in  the  block  diagram 
in  Figure  2.  The  signal -to-noi se  ratio  expression  in  Equation  (7) 
appl ies  to  r(t) . 


The  chopped  signal  will  be  approximated  by  an  amplitude-modulated 
signal  with  a  modulation  index  of  100°'  since  only  the  fundamental 
harmonic  term  of  the  signal  will  be  picked  up  by  the  audio  filter,  so 
that  the  input  to  the  square-law  device  is 


I  F 

Mi) 

SQUARE- LAW 

AUDIO 

amplifier 

■ 

DEVICE 

F  ILTER 

Receiver  system  with  square-law  detector. 


(13) 


r(t)  =  A  [l+uCOSt]cos.  t+n ( t ) 

C  a  C 

where  A  is  the  carrier  amplitude,  ..  is  the  modulation  index.  n(t)  is 

narrowband  Gaussian  noise  with  spectral  density  determined  by  the 

spectral  response  of  the  intermediate-frequency  amplifier  and  centered 

at  the  intermediate  angular  frequency  *  ,  and  j  is  the  angular  center 

C  d 

frequency  of  the  audio  filter  which  is  equal  to  the  angular  frequency 

of  the  chopped  signal.  Typical  values  for  ^  /2-r  and  _  / 2-  are  1  kHz 

a  c 

and  1  GHz,  respectively. 

The  signal -to-noise  power  ratio  at  the  output  of  the  audio  filter 
for  such  a  system  and  the  input  as  given  in  Equation  (13)  is  given  by  [3] 

( E/N)  =  „2(S./N)2  ;l  +  (S/N)(2+u2) --1  ,  (14) 

v  d 

where  (S/'N)  is  the  pre-detection  si gnal -to-noi se  ratio,  B.^  is  the 
eauivalent  rectangular  bandwidth  of  the  intermediate-frequency  ampli¬ 
fier,  and  8  is  the  bandwidth  of  the  audio  filter.  Finally,  for  a 

d 

modulation  index  of  100c',  Equation  (14)  becomes 


(E/N)  = 


(S/N)2  Ba  1 
1+3(S/N) 


with  (S/N)  as  given  in  Equation  (7). 


(15) 


This  is  the  required  expression  for  systems  employing  a  chopped- 
waveform  signal  and  square-law  detection. 


2 .  Two-way  transmission  and 
square-law  detection 

In  this  system,  the  receiver  is  identical  to  that  of  the  previous 
system.  This  time,  however,  a  target  is  used  and  the  reflected  signal 


8 


from  this  target  is  received.  Thus,  the  processed  signal -to-noise 
ratio  is  again  given  by  Equation  (15),  but  with  (S/N)  as  given  by 
Equation  (12). 


3 .  Chopped  signal  transmission  and 
quadrature  receiver 

This  system  also  consists  of  the  transmission  of  a  chopped 
signal  over  a  two-way  path  with  the  use  of  a  target.  In  contrast  to 
the  system  discussed  previously,  a  quadrature  detector  is  used  in  the 
receiver  as  shown  in  the  block  diagram  of  Figure  3,  where  u(t)  is  the 
unit  step  function,  which  is  introduced  here  only  to  show  that  the 


integrators  are  reset  at  t=0,  so  that  earlier  signals  are  disregarded. 
If  the  integrators,  which  are  discontinously  operating  devices,  are 
replaced  with  equivalent  continously  operating  filters  as  discussed 
below,  then  the  u(t)  should  be  omitted,  thus  the  signals  at  the  input 
to  the  detector  are  identical  for  Figures  2  and  3.  For  mathematical 
convenience,  the  integrators  will  be  used  in  deriving  the  output  signal 
to-noise  ratio. 


Following  Schwartz[4],  to  obtain  the  output  signal -to-noise  ratio 
first  a  single  correlation  detector  will  be  considered  as  shown  in 
Figure  4a.  The  output  of  the  integrator  is  given  by 

,T 

y(T)  =  •  r(t)s(t)dt  .  (16) 

J0 

Now,  consider  the  matched  filter  arrangement  of  Figure  4b.  The  output 
of  the  filter  at  t=T  is  given  by 

y(T)  =  '  r(t)u(t)h(T-t)dt  .  (17^ 

J  -oo 

By  the  definition  of  the  unit  step  function,  and  the  causality  require¬ 
ment  for  the  filter,  we  have 


u(t)  =  0  for  t<0 


(18) 
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Figure  3.  Quadrature  receiver. 


10 


Figure  4.  a)  Correlation  detector 
b)  Matched  filter. 


and 


h(T-t)  =  0  for  t>T  ■  (19) 

Combining  Equations  (18)  and  (19)  with  Equation  (17),  one  obtains 
rT 

y(T)  =  r(t)h(T-t)dt  .  (20) 

•'0 

If  the  filter  impulse  response  satisfies 


h(T-t)  =  s ( t )  for  0<t<T  , 

(21a) 

or  equivalently. 

h(t)  =  s(T-t)  for  t>0  , 

(21b) 

then  y(T)  is  the  same  as  in  the  correlation  detector, 
signals  in  Figure  4a  are. 

Explicitly,  the 

r(t)  =  A  cos..:  t 

(22) 

s(t)  =  C  cos^t 

(23) 

and,  from  Equation  (20),  the  filter  impulse  response  which  makes  it 
equivalent  to  the  integrator  must  be 
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(26) 


h(t)  =  C  cos[2^n-ojct]  =  C  cosu^t  =  s(t)  , 


and  imposing  the  conditions  in  Equations  (18)  and  (19),  one  obtains 


h(t)  =  C  cosu.'c t  u(t)u(T-t) 


(27) 


The  transfer  function  of  the  filter  given  by  the  Fourier  Transform  of 
its  impulse  response  is 


-T 


H(uj)  =  C  !  cosw  t  e^1  dt 
Jq  c 


(28) 


Carrying  out  the  integration,  one  obtains 


T  .  1  :^V|  ■  t  .  Jj~ 

p-r  j  (  j+w  ) sin  —x  T  j (j-u _  )o"  sin  « 

HU^e  c2  - -4 - +  e  c2  - i 


(29) 


c' . 


Since  the  system  shown  in  Figure  3  is  linear,  the  signal  and  noise 
components  of  the  output  may  be  treated  separately.  The  desired  signal 
at  the  output  is  obtained  directly  from  Figure  4a  and  Equations  (22) 
and  (23)  as 


y(  T) 


f1 

^0 


AC  cosuj  t  dt  = 
c 


ACT 

~r 


(30) 


which  leads  to  the  expression  for  the  signal  power 

■  (4 


(31) 


The  noise  power  is  given  by  [5] 


P 


n 


1_ 

2tt 


dw 


(32) 
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where  S  (w)  is  the  power  spectral  density  of  the  noise  component  of 
y(t)  which,  in  terms  of  the  power  spectral  density  of  the  noise  com¬ 
ponent  of  the  input  is  given  by 

Sny(u)  -  Snr(1JJ)lH(oj)!2  .  (33) 

Since,  for  a  filter  with  a  rectangular  frequency  response,  with  band¬ 
width  B  and  with  unit  amplitude,  the  noise  power  at  the  output  is  given 
by 


P„  =  FkT  B 
n  o 


:3*' 


then,  assuming  a  flat  input  noise  power  spectral  density,  from  Equations 
(32)  and  (34),  one  obtains 


FkT 


(35) 


For  H(w)  as  given  in  Equation  (29),  the  noise  power  at  the  output  is 
given  by 


Pn-kSnr{  :«<4l2d„ 
J  -00 

2tt 


(36) 


Assuming  w  >>  ,  and  using  Equation  (29),  Equation  (36)  becomes 

.  2(^c) 


1  C2T2 
P  =  —  S  — — 
n  w  nr  4 


si  i 


-)T 


r  r  u"uc\l . 

_  2  h 


du 


(37) 


Changing  the  variable  of  integration,  one  obtains 

sinx  y 


1  C2T2  2 
P  =  1  s  •  $ 

n  tt  nr  4  T 


-CT'  * 


dx 


(38) 


2tt 


Imposing  wc  >>  j-  once  more,  and  evaluating  the  integral  yields 


r-t1 
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p 

n 


I  sinx 

V  X 


dx 


2 

S  Ci 
nr 


(39) 


Combining  Equations  (35)  and  (39),  the  noise  power  is  given  by 


P 


n 


FkT 


C2T 


(40) 


An  input  signal  with  unknown  phase  may  be  expressed  as 


r(t)  =  A  cos(u)ct+e) 


(41) 


where  9  is  a  random  variable  with  a  uniform  probability  density  function. 
The  desired  signal  power  at  the  output  of  the  system  shown  in  Figure  3 
is  again  given  by 


P 


s 


.  2 

ACT  c 

'  2 


(42) 


but  for  the  noise  power,  this  time  one  obtains 


P  = 
n 


2FkT 


C2T 
o  ~3~ 


(43) 


Therefore,  the  signal-to-noise  ratio  at  the  output  of  this  system  is 
given  by 

A  2  1  B  ■  r 

(E/N)  =  Ps/Pn  =  r'idh^'T  =  S*ir-T  =  <S/N)TBif  (44) 
with  (S/N)  as  given  in  Equation  (12). 

This  is  the  signal-to-noise  ratio  expression  for  a  quadrature 
receiver  after  signal  processing  has  been  employed. 
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4 .  Pulsed  source  and  incoherent 
pulse  integration 

This  system  is  a  typical  range-gated  radar  utilizing  a  pulsed 
source.  It  may  be  used  as  a  transmi ssometer  by  measuring  the  reflec¬ 
tion  from  a  calibrated  corner  reflector  target.  At  the  receiver,  the 
signal  is  heterodyned  to  an  intermediate  frequency,  then  it  is  detect¬ 
ed  and  then  integrated  with  the  use  of  a  box-car  device.  A  box-car 
device  is  a  sample-  and-  hold  circuit  which  stretches  each  pulse  in 
time  so  as  to  cover  the  entire  pulse  repetition  period.  This  device 
helps  reduce  harmonics  of  the  pulse  repetition  frequency. 


Since  the  receiver  system  employs  post-detection  integration 
rather  than  pre-detection  (or  coherent)  integration,  RF  phase  infor¬ 
mation  is  lost  in  the  detection  process.  The  non-linear  nature  of  the 
detector  creates  a  rectifying  effect  and  harmonic  distortion  which 
reduces  available  signal  power.  If  n  pulses  of  the  same  signal -to- 
noise  ratio  were  integrated  by  an  ideal  pre-detection  integrator, 
after  integration  the  resulting  power  signal -to-noi se  ratio  would  be 
exactly  n  times  that  of  a  single  pulse.  If  post-detection  integration 
is  employed,  however,  a  .smaller  signal -to-noi se  ratio  would  be  observed. 
Thus,  for  the  signal -to-noi se  ratio  after  integration,  one  can  write 

(E/N)  -  n'n  ,  (45) 


where  n1  is  the  integration  efficiency  factor.  The  efficiency  factor 
in  Equation  (45)  is  not  the  same  as  the  efficiency  factor  n  as  Marcum[6] 
defines  it.  The  defining  equation  in  his  work  is 

;  s  \ 


n 


vN 


(46) 


wrrere  ( S/N ) ^  is  the  value  of  the  signal -to-noise  ratio  of  a  single 
pulse  required  to  yield  a  given  probability  of  detection,  and  ( S/N ) n 
is  the  value  of  the  signal -to-noise  ratio  of  a  single  pulse  required 
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to  achieve  the  same  probability  of  detection  when  n  such  pulses  are 
integrated.  This  factor  has  been  calculated  by  Marcum[6]  and  plotted 
by  Skol n i k[2] .  The  curves  shown  in  Figure  5  are  integration  improve¬ 
ment  factor  I  =  nn,  and  integration  loss  L  =  10  log  both  versus  n, 
the  number  of  pulses  integrated. 


Let  (S/N) n  in  Equation  (46)  be  equal  to  (S/N)  in  Equation  (45). 
Since  n  such  pulses,  when  integrated  give  the  same  probability  of 
detection  as  one  with  a  signal -to-noise  ratio  of  (S/N)-|,  this  value 
should  be  approximately  equal  to  the  combined  signal -to-noise  ratio 
of  n  pulses,  namely  (E/N)  in  Equation  (45).  Thus,  V  is  approximately 
equal  to  n.  Although  the  definition  of  efficiency  used  in  this  thesis 
and  applied  for  the  calculation  of  the  processed  signal -to-noi se  ratio 
does  not  involve  probability  of  detection,  i.e.,  it  is  assumed  to  be 
independent  of  this  factor,  the  curves  presented  here  are  still  appli¬ 
cable  to  the  calculation  of  interest  because  a  study  of  the  curves 
reveals  that  the  efficiency  factor  is  not  heavily  dependent  on  the 
threshold  level  (i.e.,  the  probability  of  detection).  Since  the 
primary  purpose  of  the  calculations  was  to  find  the  range  at  which 
(E/N)  has  a  value  of  10  dB,  the  curve  corresponding  to  ( S/N ) i = 1 0  dB 
in  Figure  5  should  be  used.  The  curve  for  which  false  alarm  number 

4 

nf  is  10  and  probability  of  detection  Pd  is  0.50  is  such  a  curve. 
Thus,  using  Equations  (45)  and  (12),  one  obtains 


(E/N)  = 


PtGtGrA2aLtLrn'n 

(4*)3R4kT  BF 
o 


1Q-2*R/10 


(47) 


where  n'  may  be  determined  using  Figure  5. 


5.  Pulsed  source  and  quadrature 
receiver 

This  system  also  is  a  radar  with  a  pulsed  source.  The  receiver 
is  a  quadrature  receiver  as  shown  in  Figure  3.  This  time  however,  the 
receiver  is  range-gated,  i.e.,  the  multiplier  input  shown  in  Figure  4a 
is  given  by 
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S(t)  =  C  cos^t  p(t) 


(43) 


where  p(t)  is  a  train  of  pulses  of  width  separation  T  .  which  are 
contained  in  the  time  interval  (0,T).  The  desired  signal  output  at 
t=T  is  given  by 

y(T)  =  1  r(t)p(t)s(t)  dt  =  .  (49) 


so  that  the  desired  signal  power  is 

p  :ACTd  2 
s  ‘  'v  2  ' 


(50) 


where  d  is  the  duty  cycle  of  the  signal  given  by 


d  =  ^  .  (51) 

r 

Again,  if  T,  and  t  are  chosen  conveniently,  the  impulse  response  of 
the  equivalent  filter  of  Figure  4b  is  given  by 


h(t)  =  C  cosu^t  p(t)  .  (52) 

It  the  filter  were  matched  to  only  one  RF  pulse  of  duration  x,  follow¬ 
ing  the  derivation  of  Equation  (40),  for  the  noise  power  at  the  output, 
one  obtains 

FkT  C2x 

PT  =  —j—  ■  (53) 


Thus,  for  the  noise  power  at  the  output  of  the  filter  with  its  impulse 
response  as  given  in  Equation  (52),  one  obtains 


P 


n 


FkT  C2rn 
o 

4 


(54) 


where  n  is  the  number  of  pulses  integrated  given  by 
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-  - 1_ 

Tr 

Using  Equations  (5C),  (51},  (54)  and  (55),  and  following  the  sare 
argument  leading  to  Equation  (44),  for  the  signal -to-noise  ratio  at 
the  output  of  the  range-gated  version  of  the  system  s'mwn  in  Figure  3, 
one  obtains 


(E/N)  =  (S/N)  T  3if  d  ,  (55) 

where  (S/N)  is  as  given  in  Equation  (12). 

C.  Attenuation  Cue  to  Rain  and  ron 

The  path  attenuation  parameter  a  appearing  in  Sections  A  and  3 
of  this  Chapter  has  a  value  of  approximately  0.3  dB/kn  [1]  near  94  go¬ 
at  sea  level  in  clear  air.  This  attenuation  is  due  primarily  to 
molecular  absorption  by  oxygen  and  water  vapor.  The  attenuation 
increases  very  markedly  during  rain  and  possibly  during  fog.  What  was 
sought  here  was  therefore  an  expression  for  attenuation  effects  on 
the  path  of  propagation  due  to  rain  or  fog.  Rain  an:  fog  are  assumed 
to  consist  of  spherical  water  droplets  of  complex  refractive  index  to 
account  for  dissipation  and  scattering  of  electromagnetic  enemy. 
Scattering  from  a  sphere  has  been  studied  by  M i e [ 7 ] ,  and  the  theory  is 
reviewed  by  Kerker[8].  Appendix  A  gives  pertinent  parts  of  the  theory 
cast  into  a  more  modern  form  with  the  use  of  vector  and  scalar  poten¬ 
tial  s. 

The  attenuation  factor  caused  by  a  density  of  smerical  droplets 
may  be  found  using  a  first  order  approximation  through  [9] 


a 


r 

4.343  i  CrY,(a)N(a)  da 
j0  -a  i 


(57) 


20 


wnere  ;  is  the  specific  attenuation  in  (dB/m),  a  is  tne  radius  of  tn- 
droplets  in  (m),  N  is  the  density  of  droplets  with  radius  a  in  -4 — 

rTW-m 

and  Q^-j.  is  the  extinction  efficiency  factor  given  by  Equation  (A59) 
in  Appendix  A. 

D .  Summary 

The  sianal -to-noise  ratio  expressions  before  signal  processing 
have  been  derived  in  Section  II. A  for  two-ended  systems  and  radar  type 
systems  as  given  in  Equations  (7)  and  (12). 

In  Section  II. 3,  expressions  for  signal-to-noise  ratio  for  five 
types  of  instrumentation  have  been  given  by  Equations  (15),  (44),  (47) 
and  (56).  The  pre-detection  signal-to-noise  ratios  (S/N).f  in  these 
equations  are  those  derived  in  Section  II. A. 

In  Section  II. C,  Equation  (57)  gives  tne  specific  attenuation  due 
to  rain  or  fog,  as  a  function  of  drop  size  distribution,  for  use  in 
the  equations  developed  in  Section  II. A. 
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SECTION  III 
RESULTS 


A.  System  Calculations 

Values  of  the  processed  signal -to-noise  ratio  were  calculated  as 
a  function  of  range  for  the  five  types  of  instrumentation.  The  range 
R  was  varied  from  500  meters  to  8  kilometers  in  100  meter  steps,  and 
the  appropriate  one  of  the  Equations  (15),  (44),  (45)  or  (56)  and  (7) 
or  (12)  in  Section  II  were  used.  In  the  calculations,  the  propagation 
path  is  assumed  to  be  uniform  or  to  consist  of  horizontal  layers  of 
homogeneous  specific  attenuation  each;  in  the  latter  case,  the  overall 
effect  of  these  layers  is  taken  into  account  by  replacing  (*R)  in  Equa¬ 
tions  (7)  or  (12)  by  £  -t  R . ,  where  *.  is  the  specific  attenuation  of 

i  1  1 

the  ith  layer  and  Ri  is  the  one-way  length  of  the  path  through  the  ith 
layer. 


Also,  routines  were  developed  to  calculate  and  plot  maximum 
attainable  range  as  a  function  of  the  specific  attenuation  a.  First, 
the  value  of  the  signal -to-noise  ratio  was  calculated  using  one  of 
Equations  (15),  (44).  (45)  or  (56)  for  a  processed  signal -to-noi se 
ratio  of  10  dB.  Then.  .  was  varied  from  30  dB/km  to  0  in  -0.1  dB/km 
steps.  Using  Equation  (7)  the  following  equation  was  obtained  and  used 


for  iterations  to  solve  for  R 

r 


Rn+1  =  r  !log10 


1 


PtGtV  Vr 

(4*)2kT  BF(S/N)  | 


i  '  2  1O910Rn 


(58) 


where  S/N  is  that  signal -to-noi se  ratio  at  the  output  of  the  inter¬ 
mediate  frequency  amplifier  which  yields  a  processed  signal -to-noi se 
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ratio  of  10  dB.  Similarly,  fcr  systems  of  the  radar  type,  by  use  of 
Equation  (12),  one  obtains 

(  r  .  .  ?  ^  ) 


n+1  2  i 


10  , 
= T- 


Wr'Vr’ 

(4  • ) 3kTQBF (S/N) J 


-  4  log,  R 
ylo  n 


For  t=30  dB/km,  the  iterations  were  started  with  R  =lm.  For  other 
values  of  i,  the  solution  corresponding  to  the  preceding  a  value 
R(  i+O.l)  was  used  for  R  The  iterations  were  stopped  when  R  ,  ~R 
lm.  This  method  was  found  to  converge  for  large  values  of  ■ .  A  differ¬ 
ent  method  was  used  for  smaller  values  of  x,  with  which  the  first  scheme 
(defined  by  Equations  (58)  and  (59))  fails.  This  was  tested  by  looking 
for  unreasonably  large  values  of  R.  For  this  case,  the  R  appearing  in 
the  denominator  in  Equation  (7)  was  used  on  the  left-hand  side,  leading 
to  the  equation  used  for  the  iterations  to  solve  for  the  range  for  the 
two-ended  system 


,  W'  Vr 

J  (4,)2kTQBF(S/N) 


-  iR  .'20 

,  n 


Similarly,  for  systems  of  the  radar  type,  Equation  (12)  leads  to 


/ - 5 - 

•  P.  6.  G  •  l  ir  - 
t  t  r  t 

V  (4-)3kT0BF(S/N) 


-2  ,R  /40 


Again,  R(  i+0.1)  was  taken  for  the  first  value  of  R  to  be  used  on  the 
right-hand  side  of  the  equation  and  the  same  test  as  for  large  ■  was 
made  to  stop  the  iterations. 


Throughout  the  calculations,  the  following  values  were  used: 
Gt  =  47  dB, 

Gr  =  53  dB 
Pt  =  200  mW 
U+L  =  11  dB, 

F  =  9  dB, 

2 

a  =  125  m  (when  applicable). 
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where  all  the  symbols  carry  the  same  meanings  defined  in  Section  II. A. 
The  gain  values  correspond  to  12"  and  24"  diameter  para¬ 
boloids.  The  transmitted  power  value  is  typical  of  an  IMPATT  diode 
source.  The  calibrated  target  is  a  corner  reflector;  this,  as  well  as 
the  system  loss  and  noise  figure  values  correspond  to  a  reasonable 
experimental  system.  These  parameters  are  inputs  read  by  the  programs 
which  would  run  equally  well  with  different  values.  They  can  be 
changed  easily. 

B.  Results  and  Discussions 

Figures  6  through  9  relate  to  the  system  with  physically  separate 
transmitter  and  receiver  and  with  a  square-law  detector  at  the  receiver 
end  (Section  II.B.l).  In  all  of  these  figures,  B  ,  the  audio 
filter  bandwidth  is  1  Hz.  Figures  6,  7  and  8  show  signal -to-noise  ratio 
versus  range  for  different  values  of  B^  with  a,  the  specific  attenua¬ 
tion,  as  a  parameter.  It  may  be  observed  from  these  curves  that  as  :i 
increases,  the  signal -to-noi se  ratio  E/N,  at  the  system  output  not  only 
decreases,  but  its  rate  of  decrease  with  respect  to  R  increases  very 
much.  The  parameter  Bif  does  riot  affect  the  E/N  values  much  for  small 
R,  but  decreasing  this  parameter  will  yield  better  E/N  values  for  large 
R.  Figure  9  shows  the  variation  of  R  as  a  function  of  a  for  an  E/N 
value  of  10  dB.  Even  under  severe  weather  conditions,  a  range  of  3  to 
4  kilometers  can  be  achieved. 

Figures  10,  11  and  12  show  the  curves  of  signal -to-noi se  ratio 
E/N  versus  range  for  a  single-ended  system  with  the  receiver  identical 
to  that  of  the  previous  system  (Section  II.B.2)  for  =  100 
KHz,  10  MHz,  1  GHz.  Figure  13  is  a  graph  of  range  versus  specific 
attenuation  for  the  same  systems  for  E/N= 1 0  dB.  Again,  the  signal-to- 
noise  ratio  is  sensitive  to  B.^  only  for  large  values  of  R.  Comparing 
Figures  10,  11  and  12  to  the  corresponding  figures  for  the  previous 
system,  one  can  observe  that  this  system  is  much  more  sensitive  to  a. 

The  performance  of  the  system  drops  very  rapidly  for  a=5  dB/km  compared 
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Figure  6.  Signal -to-noise  ratio  versus  range  for  the 
two-ended  system  with  square-law  detector. 
Intermediate  frequency  bandwidth  Bif=100 
KHz,  audio  filter  bandwidth  Ba=l  Hz. 
a:  specific  attenuation. 
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Figure  10.  Signal -to-noi se  ratio  versus  range  for  the 
single-ended  system  with  square-law  detector 
Intermediate  frequency  bandwidth  B-jf=100  KHz 
audio  filter  bandwidth  Ba=l  Hz.  a:  specific 
attenuation . 
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Figure  12.  Sfgna7-to-nofse  ratio  versus  range  for  the 
sing1 emended  system  with  square-law  detector 
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to  1=0.  Under  heavy  weather  conditons,  the  attainable  range  is  only 
1000  to  1500  meters. 

Figures  14,  15,  16  and  17  show  results  for  the  single-ended 
system  with  continuous  wave  transmission  and  with  a  quadrature  receiver 
(Section  II. B. 3).  They  also  represent  the  similar  system  with 
a  pulsed  source  (Section  II. B. 5)  if  the  duration  of  integra¬ 
tion  T  is  replaced  by  Td  where  d  is  the  duty  cycle.  Again,  the  first 
three  figures  show  E/N  versus  range  and  the  last  one  shows  range 
versus  :»  for  the  E/N  value  of  10  dB.  As  seen  in  Figure  17,  a  range 
of  approximately  2  kilometers  can  be  achieved  with  severe  weather 
conditions. 

Figures  18,  19,  20  and  21  show  E/N  versus  range,  and  range 
versus  :t  curves  for  the  single-ended  system  which  employs  post-detec¬ 
tion  integration  (Section  II. B. 4).  In  the  first  three  figures, 
the  value  of  n,  the  number  of  pulses  integrated,  are  different.  For 
all  the  curves,  the  intermediate  frequency bandwidth  is  10  MHz.  As  is 
seen  in  Figure  21,  a  range  of  approximately  1  km  can  be  achieved  with 
high  attenuation. 

Assuming  the  intermediate  frequency  bandwidth  figure  of  10  MHz 
to  be  typical  for  the  receiver  with  the  square-law  detector,  and  since 
the  value  of  1  Hz  for  the  audio  filter  bandwidth  requires  the  signal 
to  be  processed  for  approximately  1  second,  Figures  11  and  14  are  the 
curves  to  compare  the  performances  of  the  single-ended  systems  with 
quadrature  and  square-law  detectors.  It  may  be  observed  from  these 
figures  that  the  system  with  the  quadrature  receiver  is  much  superior 
especially  for  large  R  and  large  -x  alues.  This  is  expected,  because 
for  the  system  with  the  quadrature  receiver,  the  processed  signal -to- 
noise  ratio  is  proportional  to  the  signal -to-noi se  ratio  at  the  output 
of  the  i.f.  filter  as  seen  in  Equation  (44);  whereas  for  the  system  with 
a  square-law  detector,  the  final  signal -to-noi se  ratio  will  be  smaller. 
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Figure  15.  Signal -to-noi se  ratio  versus  range  for  the 
single-ended  system  with  quadrature 
receiver.  Integration  time  T=0.1  sec. 
a:  specific  attenuation. 


35 


RRNGE  (M) 

10U0  2 OOO  3000  4 000  5000 


<\J 


3=10  dB/km 


Figure  18 


3=20  dB/km  3=5  dB/km 

Signal -to-noi se  ratio  versus  range  for  the 
single-ended  system  employing  post-detection 
integration.  Intermediate  frequency  band¬ 
width  B-jf=10  MHz,  number  of  pulses  n=100000. 
a:  specific  attenuation. 


Figure  19.  Signal -to-noi se  ratio  versus  range  for  the  single 
ended  system  employing  post-detection  integration 
Intermediate  frequency  bandwidth  B-jf=10  MHz, 
number  of  pulses  n=1000.  *:  specific 

attenuation. 
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Figure  20.  Signal -to-noise  ratio  versus  range  for  the 
single-ended  system  employing  post-detection 
integration.  Intermediate  frequency  band¬ 
width  B -j f  =  1 0  MHz,  number  of  pulses  n=10. 
a:  specific  attenuation. 
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as  can  be  observed  from  Equation  (15),  especially  for  small  values  of 
unprocessed  signal -to-noise  ratio  corresponding  to  large  R  or  *. 

To  compare  the  two  radar-type  systems  with  pulsed  sources,  con¬ 
sider  first  the  curves  in  Figure  18.  The  i.f.  bandwidth  is  10  MHz 
which  implies  a  pulse  width  of  approximately  10‘7  seconds.  Assuming 

a  duty  cycle  of  0.01,  this  would  lead  to  a  pulse  repetition  period  of 
-5 

10  seconds.  Since  100,000  pulses  are  integrated,  the  entire  pro¬ 
cessing  would  take  1  second.  The  curves  for  the  single-ended  system 
with  a  quadrature  receiver  with  an  integration  duration  of  1  second 
and  a  duty  cycle  of  0.01  are  shown  in  Figure  16.  The  system  with  the 
quadrature  receiver  does  better  than  the  one  which  employs  post¬ 
detection  integration  as  is  expected  following  the  argument  given  in 
Section  II. B. 4. 

Another  observation  that  can  be  made  from  the  curves  presented 
here  is  that  changing  system  parameters  such  as  duration  of  integra¬ 
tion,  number  of  pulses,  target  cross  section  area,  etc.  prove  to  be  of 
more  value  in  attaining  larger  signal -to-noi se  ratio  for  given  range 
than  in  attaining  greater  range  for  given  signal -to-noi se  ratio.  This 
is  expected  from  Equations  (7)  and  (12).  It  can  be  observed  from  these 
equations  that  range  for  fixed  signal -to-noise  ratio  is  less  dependent 
on  these  parameters  than  is  s ignal -to-noise  ratio  for  fixed  range. 

Signal -to-noise  ratio  is  generally  a  linear  function  of  these  variables. 
Range,  however  behaves  as  a  square-root  function  for  the  two-ended 
system  and  a  fourth-root  function  for  radar  type  systems  in  clear  air 
When  attenuation  is  high,  the  exponential  terms  in  Equations  (7)  and 
(12)  become  dominant  and  the  range  is  proportional  to  the  logarithm  of 
the  parameter  in  question  for  single-ended  systems  and  half  the 
logarithm  for  the  two-ended  system. 

It  is  the  less  sensitive  dependence  of  signal -to-noise  ratio  on 
range  and  specific  attenuation  that  gives  the  two-ended  system  an 


42 


■ 


advantage  which  can  be  only  partially  offset  by  employing  coherent 
detection  with  a  single-ended  configuration. 

It  is  also  of  interest  to  see  how  a  pulsed  radar  with  pulses  of 
1  kW  peak  power  would  perform.  Figure  22  shows  range  versus  specific 
attenuation  for  a  radar  system  with  a  quadrature  receiver.  Figure  23 
is  for  a  radar  system  with  post-detection  integration.  Except  for 
transmitted  power,  all  pertinent  parameters  are  the  same  as  those 
given  on  page  23.  It  can  be  observed  that  with  the  quadrature  receiver, 
one  can  obtain  ranges  better  than  the  pulse  integration  receiver  by 
about  5C3  meters.  In  both  systems,  a  range  improvement  of  500  meters 
to  1  km  is  observed  for  a  values  of  15  to  30  dB/km  due  to  the  improve¬ 
ment  in  transmitted  power  from  200  mW  to  1  kW. 

C .  Rain  Attenuation 

Another  computer  program  developed  in  this  study  is  a  subroutine 
which  returns  the  value  of  specific  attenuation  for  a  given  rain  rate. 
The  rain  drop  size  distribution  was  assumed  to  be  exponential,  given 
by  [10] 

N(D)  =  N  e~'lD  [drops/ (m^-mm)]  ,  (62) 

where 

A  =  ar6  (mm”1 ) , 

Nq  =  8000  [drops/(m^-mm)] , 
a  =  4.1, 

6  =  -0.21 

where  r  is  the  rain  rate  in  min/Hr  and  D  is  the  droplet  diameter  in  mm. 
The  rain  water  refractive  index  was  taken  to  be  [11] 

n  =  2.2-jO. 74 
c 
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Figure  22.  Range  versus  specific  attenuation  for  E/N=10  dB 
for  the  radar  system  with  quadrature  receiver. 
Power  transmitted  Pt=1  Kw,  integration  time 
T=1  sec,  duty  cycle  d=0 . 01 . 
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Figure  23.  Range  versus  specific  attenuation  for  E/N=10  dB 
for  the  radar  system  with  post-detection  inte¬ 
gration.  Power  transmitted  Pt=1  Kw,  intermediate 
frequency  bandwidth  B-jf=10  MHz,  number  of  pulses 
n=l 00000. 
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The  equations  used  to  develop  this  program  are  Equation  (57)  in  Chapter 
II,  and  Equations  (A59),  (A63)  and  (A64)  in  Appendix  A.  The  calculated 
extinction  due  to  rain  is  shown  in  Figure  24. 

This  subroutine  was  applied  to  plot  maximum  attainable  range 
for  the  system  with  the  quadrature  receiver.  The  results  are  shown 
in  Figure  25.  The  routine  can  be  used  equally  well  with  the  other 
systems  discussed  in  Section  II. 

The  rain  drop  size  distribution  is  supplied  as  a  subprogram. 

By  changing  this  subprogram,  the  same  routine  may  be  used  to  calculate 
extinction  due  to  rain  with  other  expressions  or  empirical  data  for 
size  distribution.  This  is  useful  because  the  distribution  of  the 
smaller  drops  becomes  important  in  the  millimeter-wave  range,  and 
the  distributions  given  in  the  literature  differ  most  for  small  drop 
sizes.  Given  the  appropriate  size  distribution,  the  routine  can 
handle  attenuation  through  fog  as  well. 

D.  Conclusions 

Five  systems  considered  for  transmission  measurements  at  94  GHz 
were  studied.  One  of  these  systems  measures  transmission  directly  from 
the  transmitter  to  the  receiver.  The  other  four  systems  are  single- 
ended,  i.e.,  they  operate  like  a  radar  measuring  signals  that  are 
reflected  from  a  calibrated  target.  Each  single-ended  system  employs 
a  different  signal  processing  scheme. 

It  was  found  that  much  better  range  values  can  be  acquired  with 
the  two-ended  system  especially  under  severe  weather  conditions.  This 
is  mainly  because  with  this  system,  the  path  over  which  the  signal 
suffers  an  exponential  decay  is  only  half  of  what  it  is  with  the  other 
systems.  The  single-ended  systems  however  have  the  advantage  of  the 
ease  with  which  range  can  be  varied.  This  involves  transporting  only 
the  target.  If  for  this  reason  one  wishes  to  use  a  single-ended  system. 
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it  was  found  that  the  best  performance  is  a.nieved  with  a  quadrature 
receiver.  This  is  due  to  the  linear  nature  of  this  receiver  which 
becomes  important  especially  for  low  values  of  signal -to-noi se ' ratio 
at  the  input  of  the  detector. 


Figure  24.  Specific  attenuation  versus  rain  rate  at  94  GHz 
with  Marshall  and  Palmer  drop  size  distribution. 
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Figure  25.  Range  versus  rain  rate  for  E/N=10  dB  for  the 
single-ended  system  with  quadrature  receiver. 
T:  integration  time. 


APPENDIX  A 

SCATTERING  FROM  A  SI  MERE 


The  problem  of  scattering  from  a  sphere  was  first  solved  by 
Mie[7],  The  solution  has  been  reproduced  in  various  textbooks;  the 
treatment  given  here  follows  Kerker[8]  and  was  undertaken  primarily 
to  check  the  formulas  used  for  computation.  It  differs  from  Kerker's 
derivation  in  that,  instead  of  the  traditional  Hertz  potentials,  tne 
currently  more  popular  vector  potential  formulation  is  used. 

In  a  homogeneous,  source-free  region.  Maxwell's  Equations  for 
eJ"t  time  dependence  are  given  by 


-7*E  =  jiouH 

1 

(Al) 

7*H  =  j^eE 

(A2) 

7  •  FT  =  o 

(A3) 

7-E  =  0 

9 

(A4 ) 

where  u,  c  are  constants,  which  may  be  complex.  Equations  (A3)  and 
(A4)  suggest  that  the  fields  can  be  expresseo  in  terms  of  a  magnetic 
vector  potential  A  and  a  dual  electric  vector  potential  F.  A  and  F 
satisfy  the  equations[l 2] 

7*V*A  -  k2A  =  -ju)u7}a  ,  (A5) 


and 
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■ " ~ ' F  -  k2F  =  -j^7:f  ,  (« 

3  -f 

where  :  and  are  arbitrary  scalars.  The  field  quantities  in  terms 
of  the  vector  potentials  are  given  by 


E  =  -  ?< F  +  - -  7*7*A 


H  =  7xA  +  j~T  7*7xF 


A  =  r  A 


F  =  r  F 


(A  10) 


so  that  the  fields  can  be  expressed  as  a  superposition  of  two  parts, 
one  TM  to  r  and  one  TE  to  r,  where  r  is  the  outward  directed  unit 
vector.  A  convenient  way  to  completely  define  the  vector  potentials 
is  to  let 


,  a  . 


(All) 


r  3F 
.•  .f  _  r 

-J0)UJ  -  — 


(A12) 


Then,  from  Equations  (A5)  and  (A6),  one  obtains 


'Pr  13,.  *PA  +  1  '  F 

-2-  +  t —  ^  sin9  +  TTT - l 

^  r  sinft  V  '  r  sin  e  3$ 


~f  +  k2Pr  =  0  .  (A13) 


3r  r  sinn 


satisfied  by  both  Pr=Ar  and  Pp'Fp-  Note  that  Equation  (A13)  is 


equivalent  to 


(72  +  k2)  ^  =  0 


(A14) 


(tV)  ^  -  0 


(A16) 


Define 


A 

a  r 


(A16) 


f  Fr 


(A1 7) 


Then,  both  of  the  b  functions  can  be  shown  by  use  of  Equation  (A13) 
to  satisfy 


1  /  1  j  <  .  '--j\  1  3  , 

-  I  +  -w -  fr  isin9  ~T  j  +  “5 - 2 - 2 

r  V  3r2  r2sine  ’•  V  ^  r2sin^3:2 


+  =  0  ,  (A18) 


where  ^  or  **•/  .  The  field  quantities  in  terms  of  these  potentials 
are  given  by  [12] 


E  =  -V  k2'  <r*3) 

r  jwe  ‘ 


(A19) 


a  -1  3(r-j,f)  ,  I  32(r*d) 

tg  rsine  3$  jwer  3r30 


(A20) 


,  1  JrZl  *  Va) 


j  .  - rsin-  -r  • : 


W. 


H  •  4-  ^(r  a)  ,  4  '2(r,f) 


rsme 


J.'.  r  ,r  •  • 


and 


H  ,  .  1  iil^l  +  I  .  >V-f) 


'$  r  33  j^rsine  jrl- 
to  solve  Equation  (A18),  let 


,.(r,e,i)  =  R(r)0(9)i('i) 


(A21  ) 

(A22) 

(A23) 

(A24) 

(A25) 


This  yields  three  ordinary  differential  equations  given  by 


a  l>R(r)]  j  2  n(n+l ) 

.2  2 

dr  I  r 


rR(r)  =  0 


(A26) 


_J_d_  Mil 

sine  de  js  de 


+  n(n+l )  -  ^Vl  3(e)  =  0 

sin  e  ! 


(A27) 


and 


d-  +  m2^(cj) )  =  0 

dcb^ 


(A28) 


where  n  is  integer  and  non-negative,  and  m  can  assume  the  values  -n. 
-n+1*  — .  0,  1,  — ,  n-1 ,  n.  Equation  (A26)  has  the  two  independent 
solutions  known  as  the  Ricatti-Bessel  Functions  which  are  given  by 


JnUr)  - J  g  Jn+l/2^kr^ 


(A29) 
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and 


Nn(kr*  \/tF  "n+l/2<kr>  '  I*301 

where  J n+.-] / 2  anc^  ^n+1/2  are  half-integer  order  Bessel  and  Neumann 
Functions,  respectively.  The  solutions  to  Equation  (A27)  are  the 
Associated  Legandre  Polynomials  given  by 

0(9)  =  P^(cose)  .  ( A31 ) 

For  the  solutions  to  Equation  (A28),  one  obtains 

J>(  4> )  =  sin(m4>),  cos(m<j>)  .  (A32) 


Supposing  that  an  incident  plane  wave  propagating  along  the  2- 
axis  has  associated  with  it  the  electric  field  given  by 

-  *jk  2 


E1  =  x  e  0 


(A33) 


then,  for  the  scalar  potentials,  one  obtains 
i  a  _ 

<o 
and 


k  3  n=l  j<"‘n  HTHTTT  Jn(kor)pn(cose)cos< 


r*'f  '7?  J,  HTHTTT  3n(kor)Pn(cos6,s1n'( 

0 


(A34) 


(A35) 


The  scalar  potentials  for  the  scattered  fields  can  be  expressed  as  a 
similar  series 


rip 


sa  _  1  *  h (n-1 )  2n+l  a  A(2),b  „1d1 


7  l,  ^  '^*T7anHn '(lc0r)Pi(cose)cos*  ,  (A36) 
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’  •  Xi: 


r*Sf  '  -  ~2  J  j1"'1’  STHTTTb„  H<2>(kor)pJ(cos9>sint  ,  (A37) 

o  n=1 

while  potentials  corresponding  to  the  fields  transmitted  into  the  sphere 
are  given  by 

“  ■  ^2  j,  ■*tn_1)  HT^TT  cn  dn(klr)Pi(cose)c°s*  , 


r>P 


(A38) 


n 


tf  1  r  •  (n-1 )  2n+l  .  ,  N  ,  \  . 

"77  J,  ]  nTn+TT  dn  Jn(klr)Pn(c05-)s,n: 


(A39) 


where  kQ  and  are  the  wave  numbers  in  air  and  in  the  sphere, 
respectively.  The  Ricatti-Hankel  function  is  defined  by 


n|;2)(kr)  =  Jn(kr)  -  jNn(kr)  -  H^^kr) 


n 

,(2) 


(A40) 


where  H^+j are  the  half-integer  order  Hankel  Functions  of  the  second 

kind.  For  ipr ,  only  Jn(k1r)  have  been  chosen,  because  Nn(k^r)  are  not 

finite  at  r=0.  For  -j»s,  functions  which  satisfy  the  radiation  condition 

as  r  approaches  infinity  are  needed  and  therefore  H'  '(k  r)  have  been 

n  o 

used.  The  coefficients  a  .  b  ,  c  and  d„  can  be  found  using  the 

n  n  n  n  J 

boundary  conditions  at  r=a,  the  surface  of  the  sphere,  together  with 


Equations  (A19)  through  (A24)  yield 

^  [r(»iat»sa»  =  Ijt  [r*ta]  , 

(A41) 

^  [rf*1  Vf>]  »  ^  [r,tf]  , 

(A42) 

/  ia  sa,  .  ta 

jojeor(t|)  +iji  )  =  joie^rij) 

(A43) 

jwuQr  ( ii1  ^+'}>S^)  = 

( A44) 
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and 


(A52) 


where  k  is  the  unit  vector  in  the  direction  of  propagation  and  Z  is 
the  free-space  impedance  given  by 


Z 


c 


(A53) 


Define  the  scattering  cross-section  of  the  sphere  as 


Refd^H5*)-"  ds  ,  z  Rd 
SCA  -  Ff  3i,,  C 


(E^H5*)  *n  ds 


(A54 ) 


Re{E  ><H 


In  arriving  at  the  last  equality,  the  specific  incident  field  of 
Equation  (A33)  was  employed,  E$,  Hs  are  the  scattered  fields  corre¬ 
sponding  to  this  incident  field.  Since,  in  general,  the  refractive 
index  of  the  sphere,  nc,  may  be  complex,  some  energy  will  be  dissipated. 
Thus,  define  the  extinction  cross-section 


CEXT  "  CSCA  +  CABS  ’  (A55) 

where  C^g<.  is  the  absorption  cross-section  of  the  sphere.  Consider 
the  integral 


ZcRe  o  [(E1+ES)x(H1+Hf,)*].n  ds 


r2ir 


(eV  -E‘h:  )tc  sine  de  d4> 

u  ®  q  e 


+  Z  Re 
c 


Jo 


(-^7T  ,  c  i  *  1  C*  c  i  *  i  c*  ? 

<eV  +e!h*  -EV  -E1  )r  sine  de  d* 

g  9  9  84>  <?>  8  ?6 


+  Z  Re 
c 


IT  c2ir 


f  11  (EV*-EsH**)r2  Sine  de  d* 

gjg  9  (j)  <(>  6 


( A56 ) 


Considering  that  the  surface  on  which  the  integrals  are  being  evaluated 
does  not  enclose  any  sources,  the  first  integral  on  the  right-hand  side 
which  stands  for  the  net  outflux  of  energy  in  the  unperturbed  incident 
wave  must  be  equal  to  zero.  The  third  term  obviously  is  If 

energy  is  to  be  conserved,  the  sum  of  the  three  terms  must  be  equal  to 
-Cabs>  Since  the  sphere,  absorbing  energy,  is  the  only  energy  sink. 
Therefore,  the  second  term  must  equal  -C^y.  This  integral  has  been 
evaluated  and  yields  [7] 


CEXT-?r  „I,  <2"+')[MVb„)] 

The  corresponding  efficiency  factor  for  extinction  (dimensionless) 
is  defined  as 


^EXT  2 
?ra 


(A58) 


From  Equations  (A57)  and  (A58),  one  obtains 


VxT  =  - S  l  (2n+l  )[Re(a  +b  )] 

EXT  (k  a)2  n=l  n  n 

o 


(A59) 


The  coefficients  an  and  bn  may  be  expressed  in  terms  of  the  Spherical 
Bessel  and  Hankel  Functions  using 


=  I  ?  Jn+l/2^  '  z  Jn^ 


(A60) 


fl‘2)(Z)  -Jp  H<2>/2(z)  .  z  h<2)(z) 


(A61) 


The  derivatives  in  the  expressions  for  an  and  bn  can  also  be  eliminated 
by  use  of  the  relationship  given  by  [9] 


5T9n(z)  =  2n+l  tngn_i  (z)-(n+T  )gn+1  (z)3 


(A62) 
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where  gn(z )=jn(z )  or  h^  ' {z) . 


Using  these  Equations.  (A45)  and  (A46)  become, 


{Vkoa)><kla>  *  2n+T  kla  Jn-l(kla)  -  2STT  kla  Vl(kla>] 


a  =  - 

n  i 


<  Jn(k1a)_-Jn(koa)  +  2^1  Kp*  VlM  ~  M.  kod  WVjjj  . 

|hn2)(koa)  Jn ( kl a )  +  2rHT  kla  jn-l^kla^  "  2n?T  kla  Jn+1(k1a) 

"nc  Jn(kla)|hn2)(koa)  +  2n+T  koa  hn-!(koa)  '  2riTT  koa  hn+!(koa)'j 

(A63) 


lVkoa)  Vkla)  2n+l  kla  Jn-1(kla^  "  2n+l  kla  jn+l(kla) 

<■  ~r  •> 

;ykla)Ljn^koa)  +  2n+T~  koa  Jn-1(koa)  '  2n+T  k0a  Jn+1(koa)jj 
|hn2)(koa)^n^kla^  +  2n+I  V  jn-l(kla)  '  2n+T  kla  Jn+l(kla)j 

-JVkla)  lhn2^koa)  +  2n+T  koa  hn-l(koa)  '  2n+T  koa  hn+l(koa);  , 


(A64 ) 


APPENDIX  B 

GENERATION  OF  THE  SPHERICAL  WAVE  FUNCTIONS 


Since  the  expressions  for  Mie  coefficients  a„  and  b  involve 

n  n 

Spherical  Bessel  and  Neumann  Functions,  the  task  of  generating  these 
functions  in  a  fast  and  accurate  manner  was  encountered.  The  aim 
has  been  to  develop  a  scheme  which  could  be  used  to  evaluate  the 
functions  for  very  small  arguments  (of  order  10  )  so  as  to  enable  the 

use  of  the  routine  for  the  problem  of  modeling  attenuation  through  fog. 
It  has  been  found  that  for  arguments  as  small  as  those  of  interest 
here,  the  direct  recursion  relationship  satisfied  by  the  Bessel 
Functions  cannot  be  used  since  errors  tend  to  propagate  and  signifi¬ 
cant  digits  are  lost  after  only  a  few  recursions. 


Since  Spherical  Bessel  Functions  of  integer  order  and  Cylindrical 
Bessel  Functions  of  half-integer  order  are  related  through 


j„M  ■ 


J2x  Jn+l/2(x) 


(Bl) 


an  equivalent  problem  is  to  generate  half-integer  order  Bessel 
Functions.  For  this,  depending  on  the  magnitude  of  the  real  part  of 
the  argument,  two  methods  were  used. 


1 .  Small  Re(x) 

The  method  employed  here  makes  use  of  an  auxiliary  function  which 
is  closely  related  to  the  Bessel  Function  [13] 


A 

\) 


(X) 


r(v+D 

(x/2)v 


Jv(x) 


(B2) 
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where  r  is  the  Gamma  Function  defined  for  real  y  by 


:-(y)  =  ;  ty  e  1  dt  for  y^>0 
J0 

and  r (y )  for  y<0  may  be  found  using 


(B3) 


r(y) 


=  r(y+l 


(B4 ) 


Since  the  Bessel  Functions  satisfy  the  well-known  recursion 
relationship 


Jv+iU)  +  =  ^  JV(X) 


(B5) 


it  can  readily  be  shown  that  the  auxiliary  functions  will  satisfy 


Av(x)  Vl(x)  +  47 %+](x) 


(B6) 


It  can  also  be  shown  by  substituting  into  Equation  (B2)  that  since 


then,  the  A  functions  have  the  representation 


The  power  series  in  Equation  (B8)  has  been  found  to  converge  very  fast 
for  small  arguments.  First,  the  functions  A  (x)  for  v=61/2  and  v=59/2 
were  evaluated  employing  this  power  series  expansion.  In  the  process, 
only  four  terms  of  the  series  were  considered.  The  GMMA  subroutine  in 
the  IBM  Scientific  Subroutine  Package  was  used  to  evaluate  the  Gamma 
Functions.  Then,  using  the  recursion  formula  in  Equation  (B6)  down¬ 
ward,  Ay(x)  for  positive,  half-integer  v  less  then  59/2  were  obtained. 
Finally,  the  Spherical  Bessel  Functions  were  evaluated  by  use  of 
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Equations  (Bl)  and  (B2).  This  method  was  used  for  arguments  which 
satisfy 

Re(x)  <  3.4  . 

2 .  Large  Re(x) 

The  recursion  relationship  satisfied  by  the  Spherical  Bessel 
Functions 

Vl(x)  '  IT1  Jn(x)  '  jn-l(x)  ’  (B9> 

and  the  two  identities 

J0(x)  ,  (BIO) 


and 

j  (x)-cosx 

jl(x)  =  ~ — 7 -  (Bll) 

were  used  to  evaluate  jn(x)  for  n<nc,  where  nc  is  the  largest  integer 
smaller  than  or  equal  to  Re(x).  Then,  starting  with  an  arbitrary  value 
of  j2g(x)  and  0  f°r  J3q(x)’  the  recursion  formula  in  Equation  (B9) 
was  used  to  find  Jn(x)  for  n^nc .  Finally,  the  preliminary  values  of 
jn(x)  for  nc<n<_30  were  scaled  by  the  same  constant  so  as  to  match  the 
value  of  j  (x)  found  by  forward  recursion.  This  method  was  used  for 
the  arguments  which  satisfy  Re(x)>3.4. 

To  generate  the  Spherical  Neumann  Functions,  Bessel  Functions 
Jy(x)  for  half-integer  v<_-3/2  were  calculated  using  Equation  (B5)  and 
the  two  values  given  by 

■wx>  ■  J¥  it*-  ■  <bi2> 


and 
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J-l/2(x)  -  -M- 


cosx 

X 


(B13) 


Then,  the  Spherical  Neumann  Functions  were  evaluated  using 

n  (x)  •  (BH) 

Both  the  Sperhical  Bessel  and  Neumann  Functions  thus  generated 
using  single  precision  were  checked  against  "Tables  of  Spherical 
Bessel  Functions"[13]  and  were  found  to  be  accurate  to  five  significant 
digits. 
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